﻿ UNIVERSITATEA „POLITEHNICA” DIN TIMIŞOARA GÁL JÁNOS OLTEAN MARIUS COMUNICAŢII DE DATE ÎNDRUMĂTOR DE LABORATOR TIMIŞOARA 2004 Comunicaţii de date Lucrări de laborator Cuprins Lucrarea nr 1 Prezentarea mediului MATLAB 3 1 1 Principalele comenzi şi funcţii în MATLAB 3 1 2 Moduri de lucru în MATLAB 7 1 3 Desfăşurarea lucrării 8 1 4 Exerciţii…………………………………………………………………………………… 10 Lucrarea nr 2 Tehnici de codare digitală a semnalelor în banda de bază 11 2 1 Introducere 11 2 2 Codarea semnalelor în banda de bază 11 2 3 Prezentarea unor tehnici de codare 11 2 4 Desfăşurarea lucrării 15 2 4 1 Programul de simulare 15 2 4 2 Sarcini de îndeplinit de către studenţi 15 Lucrarea nr 3 Modelarea matematică a sistemelor de transmisiuni în banda de bază Trecerea unei secvenţe de date printr-un filtru global tip “trece-jos ideal” 16 3 1 Modelarea sistemelor de transmisiuni în banda de bază Teorema lui Nyquist 16 3 2 Simularea pe calculator a trecerii unei secvenţe de date printr-un canal echivalent de tip 18 “trece-jos ideal” 3 2 1 Metoda de echivalare a sistemelor în timp continuu cu sisteme în timp discret prin 18 invarianţa răspunsului la impuls 3 2 2 Răspunsul unui sistem în timp discret la o secvenţă de intrare x[n]………………… 18 3 2 3 Generarea cu calculatorul a unei secvenţe de date aleatoare 19 3 2 4 Răspunsul în frecvenţă al sistemului discret echivalent tip “filtru trece-jos ideal” 20 3 2 5 Secvenţa de la ieşirea sistemului numeric şi interpretarea acesteia 21 3 3 Desfăşurarea lucrării 22 Lucrarea nr 4 Transmisii de date în banda de bază fără interferenţă intersimbol, canalul fiind fără zgomot 23 4 1 Filtrul global de tip “cosinus-ridicat” Filtrul de tip “cosinus pătrat” 23 4 2 Procedura de simulare a lanţului de transmisiuni în banda de bază 25 4 3 Desfăşurarea lucrării 27 Lucrarea nr 5 Transmisii de date în banda de bază fără interferenţă intersimbol, canalul fiind cu zgomot 28 5 1 Transmisii de date prin canale cu perturbaţii 28 5 2 Procedura de simulare a lanţului de transmisiuni în banda de bază 31 5 3 Desfăşurarea lucrării 33 Comunicaţii de date Lucrări de laborator Lucrarea nr 6 Transmisii de date în banda de bază cu interferenţă controlată intersimbol 34 6 1 Transmisii de date în banda de bază cu interferenţă controlată intersimbol, canalul fiind fără zgomot 34 6 2 Transmisii de date cu interferenţă controlată intersimbol prin canale cu perturbaţii 36 6 3 Procedura de simulare a lanţului de transmisiuni în banda de bază 38 6 4 Desfăşurarea lucrării 41 Anexe 43 Bibliografie 55 Comunicaţii de date Lucrări de laborator Lucrarea nr 1 Prezentarea mediului MATLAB MATLAB (MATrix LABoratory) este un program interactiv, dezvoltat de firma Math Works Inc , fiind destinat, în special, prelucrării numerice a datelor furnizate sub formă vectorială sau matriceală MATLAB integrează calculul numeric cu vizualizarea rezultatelor şi programarea, într-un mediu flexibil şi deschis dezvoltărilor ulterioare MATLAB este destinat în mod particular prelucrării numerice (sau digitale) a semnalelor, modelării, simulării şi proiectării sistemelor numerice complexe, precum şi dezvoltării de algoritmi de prelucrare a seriilor de date, imaginilor sau câmpurilor multidimensionale de date Subsistemele de prelucrare (Toolbox-urile) extind funcţiile MATLAB-ului şi realizează operaţiuni dedicate unor domenii specifice de prelucrare numerică, cum ar fi: prelucrarea imaginilor, optimizarea, statistica, identificarea sau controlul sistemelor, studiul reţelelor neurale, a sistemelor fuzzy etc SIMULINK este un mediu interactiv pentru modelarea şi simularea sistemelor dinamice liniare, neliniare, în timp continuu, în timp discret sau a unor structuri hibride, compuse din subsisteme analogice şi digitale SIMULINK permite modelarea sistemelor dinamice cu ajutorul unor ecuaţii matematice, a unor blocuri funcţionale predefinite sau a unor blocuri definite de utilizator MATLAB este un sistem interactiv, al cărui element de bază este o matrice care nu necesită o dimensionare explicită Aceasta ajută în rezolvarea a numeroase probleme de calcul tehnic, în special pentru cele bazate pe formalismul matricial (şi vectorial), permiţând scrierea programelor într-un limbaj scalar MATLAB-ul a evoluat în ultimul timp, mai ales prin utilizarea sa în mediul universitar, ca instrument educaţional pentru cursuri de matematică, inginerie şi ştiinţă În egală măsură, MATLAB- ul se dovedeşte deosebit de eficient în cercetarea ştiinţifică şi dezvoltarea tehnologică în mediul industrial 1 1 Principalele comenzi şi funcţii în MATLAB Comenzi generale variabilele existente în spaţiul de lucru Who listează variabilele împreună cu informaţii despre ele whos listează fişierele dintr-un director cu extensiile m, mat, mex what listează dimensiunile variabilelor (de tip matrice) size precizează dimensiunile vectorilor length precizează ie on-line privind comenzile MATLAB Help documentaţ răsirea programului MATLAB exit, quit comenzi de pă Comenzi de lucru cu discul dir, chdir, comenzi similare cu cele specifice sistemului de operare DOS delete, load, save, type spaţiul de lucru pack compactează Operatori aritmetici +, -, *, /, ^ operatori aritmetici uzuali ’ operator de transpunere 3 Comunicaţii de date Lucrări de laborator Operatori relaţionali şi logici ionali uzuali , > operatori relaţ = = = operatorul de egalitate ~ = operatorul de inegalitate I logic & operatorul Ş | operatorul SAU logic ~ operatorul NU logic Comenzi generale = operatorul de atribuire ii) , separator (de argumente de funcţ i matricelor [] folosite pentru formarea vectorilor ş () folosite în expresii aritmetice : folosite în indexarea variabilelor sfârşit de linie sau suprimarea tipăririi ; indică continuarea unei comenzi pe linia următoare indică un comentariu % indică Variabile speciale şi elementare ans numele implicit al unei variabile sau al unui rezultat în virgulă mobilă eps precizia relativă pi Valoarea i, j Valoarea Inf simbolul pentru infinit NaN Not a Number iilor în virgulă mobilă flaps contor al operaţ rul argumentelor de intrare a funcţiei nargin numă rul argumentelor de ieşire a funcţiei nargout numă Funcţii matematice elementare abs valoarea absolută dăcină pătrată sqrt ră a unei variabile complexe Real partea reală a unei variabile complexe imag partea imaginară angle unghiul (faza) unei variabile complexe Conj conjugatul unei variabile complexe Sign semnul unei variabile rţiri întregi rem restul unei împă ia exponenţială exp funcţ ia logaritm natural log funcţ ia logaritm zecimal log10 funcţ 4 Comunicaţii de date Lucrări de laborator Funcţii trigonometrice ii trigonometrice uzuale sin, cos, tan, funcţ cot, sec ii trigonometrice inverse asin, acos, funcţ atan, acot, asec ii hiperbolice sinh, cosh, funcţ tanh, coth, sech ii hiperbolice inverse asinh, acosh, funcţ atanh, acoth, asech Funcţii destinate analizei datelor a unei variabile tablou Max valoarea maximă min valoarea minimă mean valoarea medie ia standard std deviaţ Sum suma elementelor unei variabile tablou a elementelor unei variabile cumsum suma cumulată tratelor elementelor unei variabile cumsqr suma pă prod produsul elementelor unei variabile tablou cumprod produsul cumulat al elementelor unei variabile Funcţii referitoare la polinoame două polinoame conv multiplică polinoame deconv împarte două rădăcinile unui polinom roots determină te un polinom din rădăcinile sale Poly construieş un polinom polyval evaluează în fractii parţiale (o fractie raţională) residue dezvoltă Funcţii referitoare la (vectori şi) matrice o matrice cu elemente nule zeros creează o matrice cu elemente unitare ones creează o matrice cu numere aleatoare distribuite uniform rand creează o matrice cu numere aleatoare distribuite normal randn creează un vector cu elemente distribuite liniar linspace creează un vector cu elemente distribuite logaritmic logspace creează valoarea unui determinant det calculează norma unui vector sau matrice norm calculează matricea inversă inv calculează 5 Comunicaţii de date Lucrări de laborator Funcţii de intrare input permite introducerea datelor de la tastatură ginput permite preluarea datelor în mod grafic, de la mouse sau cursor Comenzi grafice 1D şi 2D grafic (în planul xy) plot desenează grafic doar într o porţiune a ecranului subplot desenează un grafic cu bare bar reprezintă un grafic sub forma unei histograme hist reprezintă un grafic în coordonate polare polar reprezintă un grafic în trepte (scara) stairs reprezintă grafic o secvenţă de date (discrete) Stem reprezintă grafic în scara semilogaritmică (pe x, şi respectiv pe y) semilogx, reprezintă semilogy grafic în scara logaritmica pe x şi pe y loglog reprezintă axele de coordonate xlabel, marchează ylabel title scrie un text (explicativ) pe grafic un caroiaj pe un grafic Grid trasează o fereastră grafică figure creează close închide o figură terge graficul de pe ecran clg ş close all închide toate ferestrele grafice terge axele curente cla ş terge ecranul clr ş area suprapusă a graficelor Hold on/off permite/omite afiş (şi indică scara pe) axe axis creează Comenzi grafice 3D o suprafaţă 3D mesh desenează o suprafaţă prin curbe de nivel contour reprezintă Comenzi de control al programului ie condiţionată if comanda de execuţ else, elseif utilizate cu comanda if execuţia comenzilor for, while şi if end termină iei de un număr de ori for comanda de repetarea execuţ ionată while comanda de repetare condiţ execuţia unei bucle break termină apelul unei functii return comanda de întoarcere după teptare pause comanda de aş 6 Comunicaţii de date Lucrări de laborator 1 2 Moduri de lucru în MATLAB Lansarea în execuţie a programului MATLAB se face din WINDOWS După lansarea în execuţie, programul MATLAB intră în “modul de comanda”, afişând prompterul „ >> ” şi aşteptând introducerea unei comenzi de către utilizator Executarea unei comenzi este urmată, de obicei, fie de crearea unei variabile în spaţiul de lucru, fie de afişarea unui mesaj sau desenarea unui grafic De exemplu, comanda: >>v=0:10 va crea variabila v şi va afişa elementele acesteia pe ecran Daca nu se doreşte afişarea rezultatelor unei comenzi pe ecran, comanda trebuie urmata de caracterul „ ; ” care nu inhibă însă şi desenarea graficelor În afara modului de lucru “în linie de comanda”, în MATLAB se pot crea fişiere ce contin instrucţiuni MATLAB, numite fişiere-M (deoarece au extensia “ m”) Un program MATLAB poate fi scris sub forma fişierelor “script” sau a fişierelor “function” Ambele tipuri de fişiere sunt scrise în format ASCII, iar algoritmul care a fost implementat poate fi urmărit cu foarte mare uşurinţă, dacă se cunosc convenţiile şi sintaxa MATLAB Aceste tipuri de fişiere, obligatoriu cu extensia “ m”, permit crearea unor funcţii noi care le pot completa pe cele deja existente Prin această facilitate, MATLAB- ul poate fi extins la aplicaţii specifice utilizatorului, care are posibilitatea să scrie noi proceduri Fişiere “script” Un fişier “script” este un fişier care conţine o secvenţă de comenzi MATLAB Prin apelarea numelului fişierului, se execută secvenţa MATLAB conţinută în acesta După execuţia completă a unui fişier script, variabilele cu care acesta a operat rămân în zona de memorie a aplicaţiei Aceste fişiere nu permit integrarea în programe mari, realizate pe principiul modularizării Fişierele script sunt folosite pentru rezolvarea unor probleme care cer comenzi succesive atât de lungi, încât ar putea deveni greoaie pentru lucrul în mod iterativ, adică în modul de lucru linie de comandă Fişiere “funcţie” Dacă prima linie a fişierului-M conţine cuvântul "function", fişierul respectiv este declarat ca fişier funcţie O funcţie diferă de un fişier "script" prin faptul că poate lucra cu argumente Variabilele definite şi manipulate în interiorul fişierului funcţie sunt localizate la nivelul acestuia Prin urmare, la terminarea execuţiei unei funcţii, în memoria calculatorului nu rămân decât variabilele de ieşire ale acesteia Fişierele funcţie sunt utilizate pentru extinderea MATLAB-ului, adică pentru crearea unor noi funcţii MATLAB Forma generală a primei linii a unui fişier funcţie este: function [param ieşire] = nume funcţie (param intrare) unde: function - este cuvânt cheie care declară fişierul ca fişier funcţie (obligatoriu); nume funcţie - numele funcţiei, adică numele sub care se salvează fişierul, fără extensie Nu poate fi identic cu cel al unui fişier-M preexistent param ieşire - parametri de ieşire trebuie separaţi cu virgulă şi cuprinşi între paranteze drepte Dacă funcţia nu are parametrii de ieşire, parantezele drepte şi semnul egal nu mai au sens param intrare - parametri de intrare trebuie separaţi cu virgulă şi cuprinşi între paranteze rotunde Dacă funcţia nu are parametri de intrare, parantezele rotunde nu mai au sens 7 Comunicaţii de date Lucrări de laborator 1 3 Desfăşurarea lucrării 1 ) Se va analiza efectul execuţiei directe a unor comenzi referitoare la operaţii simple cu matrici Introducerea (directă) de la tastatură a elementelor unei matrici notată cu A se face cu: >> A = A = 2 7 9 3 1 5 8 1 2 Comanda: a = A(1:3,1) creează un vector (coloană) compus din elementele primei coloane a matricei A Răspunsul va fi: a = 2 3 8 Calculul sumei elementelor primei coloane din matricea A se face cu una din comenzile: >>sum (A(1:3,1)) sau >>sum (A(:,1)) Răspunsul va fi: ans = 13 Calculul transpusei matricei A rezultă din: >>B = A' B = 2 3 8 7 1 1 9 5 2 Extragerea elementelor diagonalei principale din matricea A rezultă din: >>diag(A) ans = 2 1 2 Analizaţi în continuare efectul următoarelor comenzi: A’, A(1,:)’, A(:), flipud(A), fliplr(A), sum(A), sum(A’), mean(A), mean(A’), sum(A,2), mean(A,2), min(A), max(A’) 2 ) Generaţi un vector C, prin întreţeserea elementelor vectorilor: A= şi B= ; 8 Comunicaţii de date Lucrări de laborator Cu secvenţa MATLAB: A= ; B= ; n=length(A)+length(B); C(1:2:n)=A; C(2:2:n)=B; se obţine rezultatul: C = 1 2 3 4 5 6 7 8 n 3 ) Se editează o funcţie genexp m, care generează o secvenţă exponenţială, x[n]=a, de lungime finită N: function x=genexp(a,n0,N) if(N =a)&(t(i) =c)&(t(i) 0  0  Considerând o succesiune de impulsuri pe linie, decalate în timp cu durata T, condiţia de interferenţă nulă între simboluri se poate scrie sub forma: (3 3) *1 ∈∀⋅= ,ℵ kkT 2F 1 din care rezultă viteza de transmitere maximă : vF == max2 T min Considerând un mesaj a(t), ce constă din N simboluri consecutive, ce apar la intervale de timp T, acesta are expresia: 1N(3 4) − () −δ= ∑ nnTta)t(a = 0n Semnalul recepţionat, r(t), în cazul unui canal ce nu introduce zgomot peste semnalul util (z(t)=0), are expresia: N−1(3 5) rt a gt nT()=⋅ () ∑ n− n= 0 fie dependent Condiţia de interferenţă nulă intersimbol presupune ca eşantionul rrkT= k(), să doar de simbolul a ş i de simbolurile ank, ≠ Notând () eşantioanele funcţiei , ki nu şnggkTk= răspuns la impuls g(t), condiţia de mai sus se verifică pentru: 17 Comunicaţii de date Lucrări de laborator ≠ 0, (3 6) gnk=∀ nk− Se observă, analizând relaţia (3 2), că în condiţiile din relaţia (3 3), se verifică întotdeauna condiţia de interferenţă nulă intersimbol din ecuaţia (3 6) 3 2 Simularea pe calculator a trecerii unei secvenţe de date printr-un canal echivalent de tip “trece-jos ideal” Pentru a putea realiza o simulare a transmisiei cu ajutorul calculatorului, plecând de la filtrul global de tip “trece-jos ideal”, trebuie să găsim un model matematic discretizat (în timp discret), echivalent cu modelul în timp continuu stabilit În acest scop vom utiliza metoda invarianţei răspunsului la impuls 3 2 1 Metoda de echivalare a sistemelor în timp continuu cu sisteme în timp discret prin invarianţa răspunsului la impuls Metoda presupune considerarea ca şi răspuns la impuls al sistemului discret echivalent secvenţa inută prin eşantionarea răspunsului la impuls g(t) al sistemului în timp gn gnT= () [] de, obţ continuu, cu un pas de eşantionare Te, ce verifică teorema eşantionării În aceste condiţii Transformata Fourier în Timp Discret, (TFTD), a secvenţei [], (), se exprimă în funcţie gnGΩ dd de răspunsul în frecvenţă () conform: Gω (3 7)   π 21   () −ω=Ω ∑ dkGG   TTT ⋅ω=Ω e   ∈Zkee Ea este o funcţie periodică în Ω, cu perioada π 2 3 2 2 Răspunsul unui sistem în timp discret la o secvenţă de intrare [] nx Notând cu [], secvenţa de la ieşirea unui sistem în timp discret, cu răspunsul la impulsul ny δ a [], la o secvenţă de intrare [], dacă sistemul este liniar şi unitate, [] ngnx n, dat de secvenţ d invariant în timp, avem: ∞ (3 8) [] [] [ ] −⋅= kngkxny ∑ −∞= k ceea ce reprezintă o sumă infinită de termeni Însă, în general, avem de-a face cu semnale de durată finită, caz în care indicele k din ecuaţia (3 8) ia valori dintr-un interval finit O cale mai rapidă de calcul a valorilor secvenţei [] de utilizarea algoritmului ny este oferită Transformatei Fourier Rapide (notat TFR sau FFT-Fast Fourier Transform), algoritm ce presupune un număr minim de operaţii necesare De fapt, aplicând TFTD, ecuaţiei (2 8) rezultă: 18 Comunicaţii de date Lucrări de laborator () () () Ω⋅Ω=Ω XGY (3 9) d ΩΩ unde () şi() reprezintă TFTD a secvenţelor []şi respectiv [] Algoritmul TFR YXnynx 2 ⋅=Ω permite calculul rapid a valorilor () cu kπ, unde M este de obicei o putere a lui 2 kYΩ Mk Calculul rapid a acestor valori se face conform: () ()() ⋅Ω=Ω (3 10) kkdkXGYΩ deci utilizând eşantioane ale spectrelor secvenţelor [] şi [] Aplicând Transformarea Fourier ngnx d Rapidă Inversă, (TFRI sau IFFT), secvenţei de valori spectrale () se obţine secvenţa răspuns, kYΩ [] ny 3 2 3 Generarea cu calculatorul a unei secvenţe de date aleatoare Avem în vedere simularea comportării unui sistem în timp continuu, caracterizat prin funcţia răspuns în frecvenţă şi răspunsul la impuls date de ecuaţia (3 2) Această analiză presupune o secvenţă de date de intrare de forma dată de ecuaţia (3 4) Ne vom rezuma doar la o transmisie binară (deci cu doar două nivele posibile ale semnalului de date), unipolară, (nivele normalizate 0 şi 1 pentru 0 şi respectiv 1 logic), şi bipolară, (nivele normalizate -1 şi 1 pentru 0 şi respectiv 1 logic) În acest scop vom găsi sistemul discret echivalent al sistemului considerat, generând şi o versiune discretă a semnalului de intrare a(t) Am văzut în paragrafl 3 1 , că echivalarea sistemelor are la bază o eşantionare cu o perioadă T a e unui semnal în timp continuu O procedură similară se aplică şi semnalelor de intrare/ieşire Astfel, () căutăm semnalul în timp discret [], semnal ce se aduce la intrare sistemului discret ednTana= echivalent Ţinând cont de structura particulară a semnalului de intrare a(t), şi anume N simboluri binare de durată T, alegem perioada T astfel încât să e nu rezulte “pierdere sau denaturare de informaţie” (scăpare de bit sau sondarea diferiţilor biţi în număr diferit de puncte) Acest lucru se asigură alegând raportul TT L= un număr întreg, de obicei o putere întreagă a lui 2 eb Considerând un număr N de biţi, tot o putere întreagă a lui 2 rezultă o lungime totală a secvenţei na în care apare o numerice, MNL=⋅ a lui 2 Se obţine o secvenţă numerică [] b, o putere întreagă d valoare de bit la intervale de timp egale cu L Dacă această valoare (0 sau 1 in cazul unipolar şi -1 b sau 1 în cazul bipolar) se generează în mod aleator, (prin ce utilizarea unor funcţii ce generează valori aleatoare), obţinem o secvenţă de date de N simboluri succesive aleatoare 19 Comunicaţii de date Lucrări de laborator Figura 3 3 O secvenţă de impulsuri unipolare Un exemplu de o astfel de secvenţă se prezintă în figura 3 3 , unde s-a considerat un număr de 4 = N=8 simboluri binare unipolare, cu 162= a secvenţei bL, rezultând M=128, (lungimea totală numerice [] na) d 3 2 4 Răspunsul în frecvenţă al sistemului discret echivalent tip “filtru trece-jos ideal” Acesta are la bază ecuaţia (3 7), rezultând o funcţie periodică, cu perioada 2⋅π în variabila Ω 2 (vezi fig 3 1a), teorema de Pulsaţia de tăiere a filtrului în timp continuu fiind ππ TF=⋅ eşantionare impune verificarea condiţiei TT1≥ Astfel, se justifică o alegere de forma e l 20, Expresia ră TT L l== spunsului în frecvenţă al sistemului discret echivalent este: eb ≥ (3 11) π  LΩ b≤ L b π ==+ 2, GGGΩ ()() ( )  dddΩΩ π π Ω ⋅+ () 001   unde α este un parametru variabil între 0 şi 1, denumit şi exces de bandă Funcţia răspuns la impuls corespunzătoare, g(t), este: (4 3) ttπ απ sin cos ∞ TTjt1ω gt G e d== ⋅ ωω ()() ∫ t2π 2π t 2α  −∞ − 1  T  T ce verifică condiţia de interferenţă nulă intersimbol: ∀≠ 0, (4 4) gn kT n k−= () () pentru α= În figura 4 2 se prezintă funcţia răspuns în frecvenţă () 05 Gω Figura 4 2 Răspunsul în frecvenţă al unui filtru de tip “cosinus-ridicat” ine funcţia de transfer global ce poartă numele de “filtru de tip cosinus pătrat”, Alegând α= 1 se obţ relaţia (4 2) devine: (4 5) ππω  2 ωω ≤  02cos = Gω ωω ()  004  ωω > 002  24 Comunicaţii de date Lucrări de laborator Funcţia răspuns la impuls are expresia: (4 6) π ttπ sin cos ∞ TTjt1ω ωω ⋅ gt G e d== ()() ∫ t22π 2π t  −∞ − 1   T   T ce verifică deasemenea condiţia de interferenţă nulă intersimbol, ecuaţia (4 4) Figura 4 3a prezintă funcţia răspuns în frecvenţă, () (ecuaţia (4 5)), în timp ce figura 4 3b Gω prezintă funcţia răspuns la impuls corespunzător, () (ecuaţia (4 6)) gt Figura 4 3a) Răspunsul în frecvenţă spunsul la impuls Figura 4 3b) Ră al filtrului “cosinus-pătrat” tor corespunză 4 2 Procedura de simulare a lanţului de transmisiuni în banda de bază Echivalând sistemul de transmisiuni în timp continuu din figura 3 1 cu un sistem în timp discret şi discretizând, deasemenea, semnalele de intrare/ieşire, se obţine structura sistemului de analizat din figura 4 4 EMIŢĂTOR CANAL RECEPTOR ad[n] Sd[n] sd’[n] rd[n] {âk} {ak Eşantionare/ ) ) GR(Ω ) MIA GE(Ω C(Ω Decizie nk Figura 4 4 Sistemul în timp discret echivalent sistemului din figura 3 1 Procedura de simulare implică: 4 2 1 Discretizarea semnalului de intrare () antionare cu pasul T Se obţine astfel o at, prin eş e l secvenţă numerică [], unde se alege TT L l==, rezultând o lungime 20,, an anT= () deeb > totală MNL=⋅ a secvenţei Alegând numărul de simboluri de transmis N ca fiind o putere b întreagă a lui 2, va rezultă o lungime totală a secvenţei M, tot o putere întreagă a lui 2 25 Comunicaţii de date Lucrări de laborator 4 2 2 Echivalarea sistemului în timp continuu cu un sistem în timp discret pe baza metodei de echivalare prin invarianţa răspunsului la impuls Se găseşte răspunsul în frecvenţă al sistemului l numeric echivalent, (), funcţie ce, în condiţiile TT L l== (ce respectă teorema GΩ 20, deb > eşantionării), va avea expresia: (4 7) LΩ 22π b = LΩΩ cos  b⋅≤0 = =+ GLGGΩ ,42π ()() ( )  dbddΩΩ  ΩΩ ω 0  T Funcţia răspuns la impuls, g(t), corespunzătoare, verifică condiţia de interferenţă nulă intersimbol: ()() ∀=− kn,Tkng≠0 (5 4) Considerând un mesaj a(t), ce constă din N simboluri consecutive, ce apar la intervale de timp T, acesta are expresia: 1N− (5 5) =− () ( ) ∑ nat a t nTδ 0n= irea acestuia se obţine Notând cu () tg ră Espunsul la impuls al filtrului de emisie, la ieş semnalul: 1N(5 6) − () −⋅= ∑ EnnTtga)t(s = 0n semnal ce urmează a fi transmis prin canal Canalul fiind idealizat, deci nu atenuează semnalul, nu-l întârzie sau dispersează, şi ţine cont doar de zgomotele ce le introduce, la intrarea filtrului de recepţie vom avea semnalul: () () () = (5 7) tztst's+ relaţie ce pune în evidenţă cele două componente ale unui semnal ce trece printr-un canal zgomotos: componenta de semnal util, s(t), şi componenta de zgomot, z(t) Dacă filtrul de recepţie este un sistem liniar atunci şi semnalul de la ieşirea acestuia va fi compus din două componente, una de semnal util şi cea de zgomot Astfel, semnalul recepţionat, r(t) se poate scrie: () () () (5 8) tztrtr+ u0= spunsul filtrului de recepţie la semnalul s(t), unde () tr reprezintă u componenta de semnal util, deci ră iar () este componenta de zgomot, deci răspunsul filtrului de recepţie la zgomotul z(t) din canal tz 0 Componenta de semnal util are expresia: 29 Comunicaţii de date Lucrări de laborator 1N(5 9) − −⋅= () ∑ nunTtga)t(r 0n = Semnalul recepţionat r(t) este eşantionat, iar pe baza eşantionului rrkT= se ia decizia k() asupra simbolului binar a transmis Semnalul recepţionat r(t) este eşantionat la momentele de timp k antioanelor obţinute (), se decide asupra simbolului transmis a tkT= rrt= kkk k, iar pe baza eş = Cum (), decizia corectă la recepţie presupune ca eşantionul de zgomot kuukakTrr= () nu depăşească un anumit prag (1/2 în cazul unei transmisii unipolare şi kTzz= k00, în modul, să 0 în cazul unei transmisii bipolare) Se defineşte raportul semnal/zgomot, RSZ, la ieşirea filtrului de recepţie la momentul t ca fiind: 0 2(5 10) () 0utr 0RSZ= P 0z irea filtrului de recepţie unde P este puterea zgomotului de la ieş 0z 2 În cazul unui zgomot alb gaussian din canal, de putere σ, puterea zgomotului la ieşirea filtrului i de recepţie este: ∞ (5 11) 2221 () ωω⋅σ=σ= dGP Riz0 ∫ 0 π 2 ∞− Semnalul util recepţionat, la momentul 0= 0t, se poate scrie: ∞ (5 12) 1 ()() () ωω⋅ω= dGSr0 Ru ∫ π 2 ∞− ω maximizeze raportul semnal/zgomot la ieşirea Se caută filtrul de recepţie optim, () RG, care să seşte ca fiind unul cu expresia răspunsului filtrului de recepţie la momentul 0= 0t Acest filtru se gă în frecvenţă: *(5 13) () () ω⋅=ω RSkG deci, filtrul de recepţie este adaptat la forma semnalului din canal Necunoscând însă, apriori, informaţia ce se va transmite prin canal, cerinţa (5 13) este dificil de realizat O soluţie mai puţin pretenţioasă, însă doar suboptimală, acceptată în practică ca şi una de compromis, este oferită de condiţia de maximizare a raportului semnal/zgomot la ieşirea filtrului de recepţie în cazul transmiterii doar a unui singur simbol de informaţie Acest lucru presupune: = a1(5 14) 0 () () ω=ω⇒ EGS ≠∀= kk,a00 30 Comunicaţii de date Lucrări de laborator ce, ţinând cont de ecuaţia (5 13), implică: 2 (5 15) ω⋅=ω⇒ω⋅=ω () () () () EERGkGGkG Condiţie satisfăcută, cu k=1, de alegerea: (5 16) ()()() ω=ω=ω GGG ER Relaţia (5 16) nu impune nici o restricţie deosebită caracteristicii de fază a celor două funcţii ωω răspuns în frecvenţă, () respectiv (), singura condiţie ce trebuie să fie verificată este acea EGRG ω de sumă nulă a celor două faze la orice frecvenţă Pentru o funcţie () reală, cum este şi cazul G nostru, putem alege: (5 17) () () () ω=ω=ω GGG RE 5 2 Procedura de simulare a lanţului de transmisiuni în banda de bază Echivalând sistemul de transmisiuni în timp continuu din figura 5 1 cu un sistem în timp discret, şi discretizând şi semnalele de intrare/ieşire, se obţine structura lanţului de transmisiuni echivalent din figură: EMIŢĂTOR CANALRECEPTOR zd[n] ad[n] S[n] sd’[n] rd[n] {âk} {ak d antionare/ Eş Ω C(Ω GR()) ) MIA GE(Ω+ Decizie nk Figura 5 2 Sistemul în timp discret echivalent sistemului din figura 5 1 Procedura de simulare implică: 5 2 1 Discretizarea semnalului de intrare (), prin eşantionare cu pasul T Se obţine astfel at e l = secvenţa [] Alegând 0,2> MNL=⋅ lLTT, rezultă an anT= () b de0= be o lungime totală n a secvenţei Alegând ,N2= va rezultă o lungime totală a secvenţei, M, tot o putere întreagă a lui 2 5 2 2 Echivalarea filtrelor în timp continuu, dimensionate conform ecuaţiei (5 17), cu filtre în timp discret Echivalarea se face pe baza metodei de echivalare prin invarianţa răspunsului la impuls Ω G ş Se găsesc funcţiile răspuns în frecvenţă ale filtrelor numerice echivalente, ()i, respectiv Ed l () Ω == ii ce, în condiţiile 02> (ce respectă teorema eşantionării), au l,LTT G, funcţ Rd0be expresia: 31 Comunicaţii de date Lucrări de laborator (5 18) Ω L bπ =≤Ω⋅ cosL 02 bΩ Ω=Ω=Ω () () () ,LGGG=  bdddRE2  π ω 0  T 2 limitat în bandă la pulsaţia π T 34 Comunicaţii de date Lucrări de laborator Răspunsul la impuls al filtrului este: 2π t(6 3) cos ∞ 41 tj ω T ()() ⋅=ωω= deGtg 2∫ ππ 42 t ∞−  − 1  T   ω În figura 6 2a se prezintă funcţia răspuns în frecvenţă, (), din ecuaţia (6 2), cu valoarea T=2, G în timp ce în figura 6 2b se prezintă funcţia răspuns la impuls, g(t), corespunzător (generat în N=128 puncte) Figura 6 2a) Răspunsul în frecvenţă Figura 6 2b) Ră spunsul în frecvenţă al filtrului “cosinus” al filtrului “cosinus” Funcţia răspuns la impuls g(t) verifică: (6 4) ()() −= 144TgTg= () ∀=⋅+− 10024,kTkTg≠ Această proprietate permite o transmisie de date cu durata 2TT= 0 a simbolurilor Astfel, dacă se consideră un semnal a(t), de intrare, ca fiind un mesaj ce constă din N simboluri consecutive, a ce apar la intervale de timp T, de forma: (6 5) {} k0 1N− (6 5) =− () ∑ 0()nat a t nTδ 0n= ω atunci, semnalul recepţionat, r(t) în cazul canalului () ce nu introduce zgomot peste semnalul C util, are expresia: 1N − (6 6) () −=nTtgatr ∑ 0)(n 0n = 0T = Eşantionând semnalul r(t) în momentele de timp kTt+, avem: 0k− 2 35 Comunicaţii de date Lucrări de laborator 1−N(6 7) TTT      000 () ⋅−−−⋅== −= + +       ∑ 10−kknkkgagaTkngar)t(r       0222=n 1−N T  0 () ⋅−−−⋅+ +=   ∑ 10−kknaaTknga   102−k,kn,n ≠= deci eşantionul prelevat la un moment t, r, depinde atât de simbolul curent cât şi de cel anterior kk transmis Într-un sistem binar monopolar de transmisiuni, simbolul a poate fi 0 sau 1, în timp ce un k sistem binar polar va emite simboluri a cu nivelele posibile -1 şi 1 Valorile posibile ale k eşantioanelor r, în cele două cazuri, se prezintă in tabelul 6 1 k Tabelul 6 1 Valorile posibile ale eşantioanelor r, în funcţie de simbolurile transmise k Sistem unipolar Sistem bipolar a a r a a r k1−kkk1−kk 0 0 0 -1 -1 -2 +1 0 +1 +1 -1 0 0 +1 +1 -1 +1 0 +1 +1 +2 +1 +1 +2 La recepţie, decizia asupra simbolului binar transmis, a, notată aˆ, se poate face pe baza kk eşantionului curent r şi a deciziei anterioare aˆ Metoda de detecţie utilizată permite deasemenea k k1− detectarea unor decizii incorecte Astfel, de exemplu, în cazul unei transmisii bipolare, o valoare r k de -2 nu poate fi urmată de o valoare 2+ facilitate însă nu elimină posibilitatea luării 1= kr Această + unor decizii incorecte, mai mult, din cauza condiţionării deciziei curente de decizia anterioară, apare posibilitatea “propagării” unei decizii incorecte (luări de decizii incorecte până la apariţia unei noi erori de decizie sau apariţia unei violări de semn) 6 2 Transmisii de date cu interferenţă controlată intersimbol prin canale cu perturbaţii Modelarea unui sistem de transmisiuni printr-o caracteristică globală de tip “cosinus”, cu idealizarea canalului permitea o analiză a metodei de transmitere cu interferenţă controlată intersimbol În realitate însă, trebuie să ţinem cont şi de influenţa canalului de comunicaţii asupra semnalului transmis prin el (atenuare, întârziere, distorsionare, adăugare de zgomot) EMIŢĂTOR CANAL RECEPTOR z(t) a(t) s(t) s’(t) r(t) {âk} {ak Eşantionare/ G(ω MIAE) GR(ω C(ω) Decizie ) + tk Figura 6 3 Schema de principiu a unui sistem de transmisiuni în banda de bază În figura 6 3 se prezintă structura unui sistem de transmisiuni în banda de bază, ce modelează un canal zgomotos În cele ce urmează vom neglija efectele de întârziere şi de distorsionare ale canalului, ţinând cont doar de zgomotul introdus de el, pe care-l vom considera ca fiind unul aditiv, 36 Comunicaţii de date Lucrări de laborator alb şi cu distribuţie gaussiană, cu dispersia constantă σ Astfel, vom considera ()1=ωC Influenţa i zgomotelor introduse de canal se pune în evidenţă prin însumarea semnalului de la ieşirea canal, () ts, cu un “zgomot”, z(t) Notând cu () răspunsul la impuls al filtrului de emisie, la ieşirea acestuia se obţine tg E semnalul: 1N − (6 8) () −=)(nTtgats ∑ 0En 0n = semnal ce este transmis prin canal Canalul fiind idealizat, deci fiind unul ce nu atenuează şi nici nu întărzie sau dispersează, la intrarea filtrului de recepţie vom avea semnalul: () () () = (6 9) tztst's+ relaţie ce pune în evidenţă cele două componente ale unui semnal ce trece printr-un canal zgomotos: componenta de semnal util, s(t), şi componenta de zgomot, z(t) Filtrul de recepţie fiind un sistem liniar rezultă că şi semnalul de la ieşirea acestuia se poate descompune în două componente, una de semnal util şi cea de zgomot Astfel, semnalul recepţionat, r(t), se poate scrie: () () () (6 10) tztrtr+ u0= unde () reprezintă componenta de semnal util, deci răspunsul filtrului de recepţie la semnalul s(t), tr u iar () ie la zgomotul z(t) din canal tz este componenta de zgomot, deci ră 0spunsul filtrului de recepţ 0T , avem: Eşantionând semnalul r(t) în momentele de timp kTt+ = 0k− 2 () () +==+= kukkkkukzrrtztrtr= 0)((6 11) N−1 TTT    000 () ⋅−−−⋅= −=+ +    ∑ kkkngagazTknga+ 10 −    n0222 = N−1 T  0 () ⋅−−−⋅+ +=+  ∑ kkkknzaazTknga+10 −  kknn1,,02 −≠= La recepţie, decizia asupra simbolului binar transmis, a, notat aˆ, se poate face pe baza kk eşantionului curent r şi deciziei anterioare aˆ Această condiţionare a deciziei curente de decizia k1−k anterioară face posibilă “propagarea” unei decizii incorecte (ce se manifestă prin luări de decizii incorecte până la apariţia unei erori noi) În practică această posibilitate se elimină prin precodarea secvenţei de date ce urmează a fi transmis Astfel, la emisie, în loc de secvenţa {} ase transmite k secvenţa {} de: kbdată (6 12) ⊕= 1− kkkbab Se defineşte raportul semnal/zgomot, RSZ, la ieşirea filtrului de recepţie la momentul t ca fiind: 0 37 Comunicaţii de date Lucrări de laborator 2(6 13) () 0utr 0RSZ= P 0z unde P este puterea zgomotului de la ieşirea filtrului de recepţie 0z ω Se caută filtrul de recepţie optim, (), care să maximizeze raportul semnal/zgomot la ieşirea RG filtrului de recepţie la momentul 0= seşte ca fiind unul cu expresia răspunsului 0t Acest filtru se gă în frecvenţă: *(6 14) () () ω⋅=ω RSkG deci, filtrul adaptat la forma semnalului din canal Necunoscând însă, apriori, informaţia ce se va transmite prin canal, cerinţa (6 14) este dificil de realizat O soluţie mai puţin pretenţioasă, însă doar suboptimală, acceptată în practică ca şi una de compromis, este oferită de condiţia de maximizare a raportului semnal/zgomot la ieşirea filtrului de recepţie în cazul transmiterii doar a unui singur simbol de informaţie Acest lucru conduce la o expresie forma: 2 (6 15) () () () () ω⋅=ω⇒ω⋅=ω EERGkGGkG În cazul în care k=1, rezultă: (6 16) ()()() ω=ω=ω GGG ER Relaţia (6 16) nu impune nici o restricţie deosebită caracteristicii de fază a celor două funcţii ωω ie ce trebuie să fie verificată este răspuns în frecvenţă, () EG respectiv ()RG, singura condiţ G reală ω aceea de sumă nulă a celor două faze la orice frecvenţă Pentru o funcţie () , cum este şi în cazul nostru, putem alege: (6 17) () () () ω=ω=ω GGG RE 6 3 Procedura de simulare a lanţului de transmisiuni în banda de bază Echivalând sistemele de transmisiuni în timp continuu din figurile 6 1 şi 6 3 cu un sistem în timp discret şi discretizând, deasemenea, semnalele de intrare/ieşire, se obţine structura sistemului de analizat din figură: TOR CANAL RECEPTOR EMIŢĂ ad[n] Sd[n] sd’[n] rd[n] E{ak{âk} şantionare/ C(Ω) ) GR(Ω MIA GE(Ω ) Decizie nk Figura 6 4 Sistemul în timp discret echivalent sistemului din figura 6 1 38 Comunicaţii de date Lucrări de laborator şi respectiv: EMIŢĂTOR CANAL RECEPTOR zd[n] ad[n] S[n] sd’[n] rd[n] {âk} {ak d Eşantionare/ ) MIA G(Ω ) + GR(Ω E) C(Ω Decizie nk Figura 6 5 Sistemul în timp discret echivalent sistemului din figura 6 3 Procedura de simulare în cazul canalului fără zgomot implică: 6 3 1 Discretizarea semnalului de intrare (), prin eşantionare cu pasul T Se obţine astfel o at e l = secvenţă numerică [], unde se alege ,l,LTT02> rezultând o lungime an anT= () debe0= totală MNL=⋅ a secvenţei Alegând numărul de simboluri de transmis N ca fiind o putere b întreagă a lui 2, va rezulta o lungime totală a secvenţei, M, tot o putere întreagă a lui 2 6 3 2 Echivalarea sistemului în timp continuu cu un sistem în timp discret pe baza metodei de echivalare prin invarianţa răspunsului la impuls Se găseşte funcţia răspuns în frecvenţă a sistemului l = numeric echivalent, (), funcţie ce, în condiţiile 02> (ce respectă teorema l,LTT GΩ be d0= eşantionării), va avea expresia: (6 18) Ω π  bL Ω=≤Ω⋅ 2cosL  0b ()() ( ) π+Ω=Ω =Ω 22GG,LG  ddbd  π , rezultă o lungime totală an anT= lLTT () de0= be n ei Alegând ,N2= va rezulta o lungime totală a secvenţei, M, tot o putere MNL=⋅ b a secvenţ întreagă a lui 2 6 3 7 Echivalarea filtrelor în timp continuu, dimensionate conform ecuaţiei (6 17), cu filtre în timp discret Echivalarea se face pe baza metodei de echivalare prin invarianţa răspunsului la impuls Ω Se găsesc funcţiile răspuns în frecvenţă ale filtrelor numerice echivalente, () şi, respectiv G Ed l () Ω == ii ce, în condiţiile 02>, (ce respectă teorema eşantionării), au l,LTT G, funcţ Rd0be expresia: (6 21) Ω L bπ =≤Ω⋅ cosL 02 bΩ Ω=Ω=Ω () () () ,LGGG=  bdddRE2  π 0 5, y(i)=1; else y(i)=0; end end x receptionat=[]; for i=1:N, x receptionat=[x receptionat y(i) zeros(1,31)]; end %daca peste semnalul util nu s-ar fi suprapus zgomot Y1fzg=X1 *Gr; y1fzg=ifft(Y1fzg); Pzgout=sum(abs(y1fzg-y1) ^2); Psout=sum(abs(y1fzg) ^2); RSZout=Psout/Pzgout; IRSZ=RSZout/RSZin; figure(1); plot(x); title('Simboluri de date emise si supraesantionate');hold on; stem((1:Lb:M),zeros(1,M/Lb),'r');hold off; figure(2); subplot(211); plot(frecv,Ge(1:2*(M*Omega0/(2*pi))));grid; title('Caracteristica de frecventa a filtrului radacina din cosinus patrat(Ge)'); subplot(212); plot((1:256),ge);grid; title('Raspunsul la impuls al filtrului de tip radacina din cosinus patrat'); 48 Comunicaţii de date Lucrări de laborator figure(3); plot(real(x1)); title('Semnalul util care se transmite in canal'); grid; figure(4); plot(zgomot); title('Zgomotul suprapus peste semnalul util'); grid; figure(5); plot(real(x2)); title('Semnalul perturbat de zgomot adus la intrarea filtrului de la receptie'); grid; disp('Puterea medie a semnalului este :');disp(Psemnal/M); disp('Puterea medie a zgomotului este :');disp(Pzgomot/M); disp('RSZ la intrarea filtrului adaptat este :');disp(RSZin); disp('Apasati o tasta pentru a continua '); pause; figure(6); plot(real(y1));hold on; stem((1:Lb:M),real(y1(1:Lb:M)),'r');hold off;grid; title('Semnalul la iesirea filtrului, cu sublinierea momentului de esantionare');hold on; prag=(0 5) *ones(1,M); plot( ,prag,'b ');hold off; figure(7); plot(real(y1));hold on;plot(real(y1fzg),'g-');hold on; title('O comparatie intre forma semnalului receptionat cu/fara zgomot aditiv'); legend(' : semnal afectat de zgomot',' : semnal neafectat de zgomot'); stem((1:Lb:M),zeros(1,N),'r');hold off;grid; disp('Puterea semnalului util la iesire este :');disp(Psout/M); disp('Puterea zgomotului la iesire este :');disp(Pzgout/M); disp('RSZ la iesirea filtrului adaptat este :');disp(RSZout); disp('Raportul semnal zgomot s-a imbunatatit de : ');disp(IRSZ); disp('Apasati o tasta pentru a continua '); pause; figure(8); subplot(211); plot(x); title('Simboluri emise');hold on; subplot(212); plot(x receptionat); title('Simboluri detectate'); 49 Comunicaţii de date Lucrări de laborator Lucrarea 6 Transmisiune de date in banda de baza cu interferenta controlata inter-simbol Trecerea unei secvente de date printr-un filtru de tip cosinus, canalul fiind fara zgomot clear; t=0:1/32:8-1/32; M=length(t); %parametrii programului Te=1/32; T=1; Lb=T/Te; Omega0=pi/Lb; N=M/Lb; %generarea unei secvente aleatoare unuipolare de N=8 simboluri si supraesantionarea ei cu Lb=32 x=[]; for i=1:N, sursa(i)=randn; if sursa(i) =0, decizie=1; y(1)=decizie; else decizie=-1; y(1)=decizie; end for i=2:N, y1(i)=real(x1(-Lb/2+i*Lb)); if (y1(i)-decizie)>=0, y(i)=1; else y(i)=-1; end decizie=y(i); end x receptionat=[zeros(1,15)]; for i=1:N, x receptionat=[x receptionat y(i) zeros(1,31)]; end 51 Comunicaţii de date Lucrări de laborator x receptionat=[zeros(1,15) x receptionat(240) zeros(1,31) x receptionat(16) zeros(1,31) x receptionat(16+32) zeros(1,31) x receptionat(16+2*32) zeros(1,31) x receptionat(16+3*32) zeros(1,31) x receptionat(16+4*32)zeros(1,31) x receptionat(16+5*32) zeros(1,31) x receptionat(16+6*32) zeros(1,16)] figure(1); plot(x); title('Simboluri de date emise si supraesantionate');hold on; stem((1:Lb:M),zeros(1,M/Lb),'r');hold off; figure(2); subplot(211); plot(omega,[G(129:256) G(1:128)]);grid; title('Caracteristica de frecventa a filtrului cosinus'); xlabel('omega');axis([-pi/2 pi/2 0 60]); subplot(212); plot((1:256),g);grid; title('Raspunsul la impuls al filtrului de tip cosinus '); culoare=['g','r','c','m','y','k','g','r','c']; figure(3); for i=1:N, plot( ,out{i},culoare(i));hold on; end plot( ,real(x1),'-o'); hold off; grid; axis( ); title('Raspunsul filtrului la secventa de date de la intrare'); figure(4); plot(real(x1));hold on; stem((Lb/2:Lb:M-Lb/2),y1,'r');hold off;grid; title('Semnalul la iesirea filtrului, cu sublinierea momentului de esantionare'); figure(5); subplot(211); plot(x); title('simboluri emise'); subplot(212); plot(x receptionat); title('Simboluri detectate'); %daca se utilizeaza o pereche de filtre clc; disp('In continuare vom retransmite secventa de date, utilizand o pereche de filtre de tip "square root cosine" '); pause; figure(6); plot(real(ge)); grid; 52 Comunicaţii de date Lucrări de laborator title('raspunsul la impuls al filtrului de tip "square root cosine’); figure(7); plot(omega,[Ge(129:256) Ge(1:128)]);grid; title('Caracteristica de frecventa a filtrului "square root cosine"'); xlabel('omega');axis([-pi/4 pi/4 0 10]); %calculul secventei de date dupa trecerea prin filtrul de emisie Xe=X *Ge; xe=ifft(Xe,256); figure(8); plot(real(xe)); title('Simbolurile de date dupa trecerea prin filtrul de emisie'); grid; %peste semnalul util se adauga zgomot xr1=xe+zgomot; figure(9); plot(real(xr1));hold on; grid; title('Semnalul afectat de zgomot care va fi prelucrat de filtrul de receptie'); plot(real(xe),'r ');hold off; %calculul secventei de date de la iesirea filtrului de receptie; Xr1=fft(xr1); Xr=Xr1 *Gr; xr=ifft(Xr); figure(10); plot(real(xr)); hold on; moment esant= ; stem(moment esant,real(xr(moment esant)),'r'); hold off; title('Semnalul de la iesirea filtrului cu sublinierea esantioanelor ce vor fi prelevate'); grid; %luarea deciziei asupra simbolurilor emise yr(1)=real(xr(Lb/2)); if yr(1)>=0, decizie=1; yr(1)=decizie; else decizie=-1; yr(1)=decizie; end for i=2:N, yr(i)=real(xr(-Lb/2+i*Lb)); if (yr(i)-decizie)>=0, yr(i)=1; else yr(i)=-1; 53 Comunicaţii de date Lucrări de laborator end decizie=yr(i); end %rearanjarea datelor yreceptionat=[zeros(1,15) yr(N) zeros(1,Lb) yr(1) zeros(1,Lb) yr(2) zeros(1,Lb) yr(3) zeros(1,Lb) yr(4) zeros(1,Lb) yr(5) zeros(1,Lb) yr(6) zeros(1,Lb) yr(7) zeros(1,Lb/2)]; figure(11); subplot(211); plot(x); title('simboluri emise'); subplot(212); plot(yreceptionat); title('Simboluri detectate'); 54 